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We deal with the presence of topological defects in models for two real scalar fields. We comment 
on defects hosting topological defects, and we search for explicit defect solutions using the trial 
orbit method. As we know, under certain circunstances the second order equations of motion can 
be solved by solutions of first order differential equations. In this case we show that the trial orbit 
method can be used very efficiently to obtain explicit solutions. 



PACS numbers: ll.10.Lm, 11.27.+d, 98.80.Cq 
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In the seventies there appeared a great deal of investi- 
gations concerning the presence of topological defects and 
their role in high energy physics. Particularly interesting 
issues appeared in the search for defects in models involv- 
ing real scalar fields, as for instance in the investigations 
of Refs. [[tj H D which considered models described by 
two real scalar fields. The interest has been renewed 
in the nineties, where there has appeared investigations 
dealing with systems of two scalar field having distinct 
motivations, as for instance in the case of defects spring- 
ing in the form of junctions of defects 
defects having internal structure 
investigations concern the presence of doman walls in su- 
pergravity |TT[ , in scenarios for localization of gravity 




and also as 
Other lines of 



on domain walls 12 , in supersymmetric gluodynamics, 
where nonperturbative effects may give rise to gluino con- 
densates [[L3| [uj, and also in string theory, since there 
are models in field theory which correctly describe the 
low energy world volume dynamics of branes in string 
theory ||[T§]l7j. 

A central issue in the investigation of defects in sys- 
tems involving two real scalar fields concerns integrabil- 
ity of the equations of motion. From the mathematical 
point of view the problem is hard, because one starts 
with two coupled second order nonlinear ordinary dif- 
ferential equations - see Refs. [jf], |[ [| fo r more details. 
However, in the investigations in Ref. Tl8| one proposes a 
new route, in which the mathematical barrier is simplified 
if one considers a specific class of systems. In this class 
of systems the equations of motion can be reduced to 
first order differential equations, which allows obtaining 
Bogomol'nyi-Prasad-Sommerfield (BPS) states (l^, p0[ , 
which are stable configurations that minimize the energy 
of the topological solutions. 

More recently, in Refs. 0, ^| one gets two new re- 
sults, one showing that under certain conditions pl| the 
second order equations of motion can be reduced to a 
family of first order equations, in this case ensuring that 
all the topological solutions arc BPS states, and the other 
showing that sometimes |2^] it is possible to find an inte- 
grating factor for the first order equations which allows 
obtaining all the BPS states of the system, thus unveiling 
the moduli space of topological kinks. 

As one knows, models described by two real scalar 



fields are generically described by a Lagrange density, 
which contains the usual kinetic and gradient contribu- 
tions, and is otherwise specified by a potential, in gen- 
eral a smooth function of the two scalar fields. If the 
potential has several minima, one may find defects of 
at least two different types: topological defects, which 
starts in a given minimum of the potential and ends in 
another, distinct minimum; nontopological defects, which 
starts and ends in the very same minimum of the poten- 
tial. These defects constitute orbits in the plane (4>, %), 
an orbit being a function of the two fields which can be 
represented by F(tfi,x) = 0. In Ref. H the author has 
proposed a method to help searching for explicit solu- 
tions to the equations of motion. The procedure has been 
named trial orbit method, since it relies on trying an or- 
bit F((f>,x) = 0, and checking its validity. Although the 
method may help finding explicit solutions, it is not that 
strong since it is directly related to solving the equations 
of motion, which are second order differential equations. 
In the present work, however, we revisit the trial orbit 
method, adapting the methodology to the investigation 
of BPS states. 

This possibility appears very naturally when one con- 
siders a specific class of models described by two real 
scalar fields. We follow recent work Q and references 
therein to introduce models defined by 



v^,x) = \wl + \w 2 x 



(1) 



Here W — W(<fi, \) is a smooth function of the two fields, 
and W$ — dW/d(f> and so forth. In this case the energy 
density for static field configurations (j)( x ) an d x( x ) can 
be written in the form 



£(x) 



dW 

dx 



1 (d<t> 
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1 (dx 
2\dx 



(2) 



and for the static field configurations that obey the first 
order differential equations 



dx 



dx 
dx 



(3) 



the energy is minimized to the bound Ebps — \AW\, 
where AW = W[(p(oo), x(oo)] - W[(j)(-oo), x(-oo)] 
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These solutions are BPS states, which solve the first or- 
der differential equations and the equations of motion of 
the model. 

To search for defect solutions we first investigate the 
minima of the potential, which are critical points of W 
obtained by requiring that W$ = and W x = 0. We 
suppose the model has the discrete set of n minima, 
(0dXi) 5 (02 ,X 2 ))•••) (^„>XJ- In principle, each pair of 
minima forms a topological sector, so we use the sub- 
scripts ij to identify the topological sectors. Some of the 
topological sectors may have solutions that obey the first 
order equations (||). In this case they are named BPS 
sectors, and we can identify the BPS sectors with the 
energy of the topological sector, with Eij 



AWij], for 



We eliminate the spatial coordinate in the first order 
equations (||) to write W x d(f)—W<p dx = 0. This equation 
is exactly solvable for W harmonic, for + W xx — 0. 
In this case the solution is given by the orbit F((f>, \) = 0, 
which is obtained via dF/d<j) — W x and dF/d\ = —Wq. 
In Ref. pl| one has obtained that when W is harmonic 
all the topological solutions are of the BPS type, that 
is, are solutions of the first order Bogomol'nyi equations. 
We use this result together with the above result to in- 
troduce a new result, which ensures that the equations of 
motion for topological solutions are exactly solved when 
the potential is of the form (|l|), with W harmonic. 

If W is not harmonic, we can yet search for some in- 
tegrating factor, / = I((f),x), such that the solution is 
now given by F(<j),x) — 0, with dF/dcj} = I(4>,x) W x 
and dF/dx = — !(<!>, x) W<£- The problem is that it is 
not always ease to find an integrating factor for the first 
order equation. This difficulty has led us to revisit the 
trial orbit method, to see if it can be of some help in the 
process of finding solutions to the first order equations. 
As we will show, the trial orbit method is of good use 
for searching for topological defects in the BPS sectors 
that appear from the first order Bogomol'nyi equations. 
Better than that, the trial orbit method is more efficient 
than it appears in the original work ||, since here it re- 
lies on searching for solutions of first order equations, 
circumventing the weakness that appears in Ref. ||, in 
which one deals with the equations of motion, which are 
second order differential equations. 

The trial orbit method may be seen as a procedure 
based on the three steps: 

Step 1. We sellect the BPS sector. We do that by 
supposing the pair of minima (4> i , \ i ) and (cj> , \. ) is such 
that W{$> i ,Xi) 7^ W / (<ApXj); implying that one is deal- 
ing with a BPS sector. This means that the first order 
equations has topological solutions connecting the points 

(0»>X») and O^Xj); 

Step 2. We choose an orbit. We do that by writing 
the equation of the orbit, say Fij{4>, x) = 0> an d checking 
compatibility between the chosen orbit and the minima 
that specify the BPS sector, that is, validating the state- 



ments F lj (4>.,x i ) 



and Fi, 



0: 



and the first order equations. We do that by differ- 
entiating the trial orbit. We get (dFij/d(j))(d(j)/dx) + 
(dFij I dx)(dx/ dx) = 0. We use the first order equations 
(||) to obtain 



BR, 



dFu 
dx 



(4) 



This new statement is similar to the trial orbit itself. 
Thus, we have to check compatible between the chosen 
orbit and the new statement We say that Fij((j), x) = 
is a good orbit when every pair (<kx) that solves the 
orbit also solves the new statement (Q). We then use the 
good orbit to solve the first order equations. 

We examplify the procedure with some specific inves- 
tigations. For simplicity we deal with natural units, and 
with dimensionless fields and coordinates. The first ex- 
ample is given by 



f 



w(<t>,x) = 4>--z<r-r<i>x 



(5) 



where the parameter r is real. This is the model first 
investigated in Ref. |l8|| . We use Eq. (||) to write the 
first order equations 



dx 

da 

dx 



4> ~rx' 



— — = — 2r(f>x 



(G) 
(7) 



Step 3. We test compatibility between the orbit 



We first consider r real and negative. The potential has 
minima at v\ 2 — (ilj 0), which define a topological sec- 
tor of the BPS type. An explicit solution describes a 
straight line orbit connecting the two minima; it corre- 
sponds to <p(x) = tanh(a;) and x — 0- F° r r = ~ 1 wc 
see that W is harmonic, and the model can be solved ex- 
actly. There is no orbit connecting the minima (±1,0), 
but the single straight line segment with x = 0. This 
result is valid not only for the first order equations, but 
also for the equations of motion, that is, the system has 
no topological solutions connecting the minima (±1, 0) 
for r = — 1, unless the kinks with (f>(x) = ±tanh(a;) and 
X = 0. 

We now consider r real and positive. In this case one 
can find an integrating factor and exactly solve the model 
p2[ . However, since we want to show how the trial or- 
bit method works we proceed using the methodology ex- 
plained above. The minima of the potential are now at 
fi,2(±l,0) and ^4 = (0, ±y/TJr). The model has one 
BPS sector with energy E12 — 4/3, and four BPS sec- 
tors with energy Ei 3 — En = E23 = E 2 i = 2/3. To 
find non trivial solutions we use the above procedure. 
We follow the first step to sellect the BPS sector de- 
fined by the minima (±1,0), with energy E12 = 4/3. 
The second step leads us to choose the orbit. We try 
1 + a<j) 2 + bx 2 =0. The ending points (±1,0) demand 
that 1 + a = 0, so we can write 1 — <fi 2 + bx 2 = 0. We 
follow the third step and we differentiate this orbit to 
get —<j){d<t>/dx) + bx(dx/dx) = 0. We use the first order 
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equations to get 1 — <fi 2 + r(2b— l)\ 2 = 0. This new state- 
ment is fully compatible with the orbit for b = r(2b — 1), 
that is, for b = r/(2r — 1). This finally gives the good 
orbit 



-KW 2 )(1 



s4> 2 ) X 2 + +2r 2 (j) 2 x 2 (13) 



1 - 2r 



x 2 = i 



(8) 



which connects the minima (±1,0) for < r < 1/2. 

We now use the orbit (||) to find explicit solutions. It 
allows rewriting Eq. (|^) in the form d<j>/dx = 2r (1 — <p 2 ), 
which is solved by <p(x) = tanh(2rx). The other field is 
then given by 



X(x) = ± 



cosh(2ra) 



(9) 



We notice that the limit r — > 1/2 leads to x = and 
4>{x) — tanh(x), which also solve the first order equa- 
tions (pj). This pair corresponds to a straight line segment 
connecting the minima (±1,0). It is different from the 
other pair, which connects the same minima by an ellip- 
tic segment. The presence of the elliptic orbit manifests 
the possibility of domain walls having internal structure. 
This possibility emerges after examining the masses of 
the two fields, which are: m 2 ^ — 4, to 2 = Ar 2 at the 

minima (±1,0), and rrA — Ar, to 2 — Ar at (0, ±^/l/r). 
Also, the energies of the topological defects are given by: 
for x — we have E\i = 4/3 and for <j> = we have 
E34 = i/3^/r. We see that for r € (0, 1) energy consid- 
erations favor the BPS defect to be the host defect, and 
for r > 1 the non BPS defect becomes the host defect. In 
both cases the elementary <j> mesons prefer to live inside 
the host defect, while the x mesons prefer to live outside. 

In the other BPS sectors we can consider orbits like 
r x 2 = 1 ± 4>i which requires r — 1/4. For r x 2 — 1 — </> we 
have 2rx{dx/dx) = —dcjy/dx. In this case we use Eqs. (||) 
i> — <fi 2 , which is solved by cj>(x) = 
The other field is given by x( x ) — 
For rx 2 = 1 + 4> the investigation 
>(x) = -(1/2)[1 - tanh(x/2)] and 



to obtain d(f>/dx 
(l/ 2)[l + tanh(a?/2)] . 
± v / 2[l -tanh(ir/2)]. 



is similar. We get 
X(x) = ± x /2[l±tanh(x/2)]. 

In order to further illustrate our procedure we consider 
other models, and we explore the most interesting BPS 
sectors they have. The first model is defined by 

W(0, X ) = ^-i(l + s)0 3 + ^ 5 - r(f>x 2 (10) 
3 5 

It can be seen as an extension of the former model, which 
includes the fifth order power on <fi. We suppose the 
parameters r and s are real. The first order equations 
are 

^ = l-(l + ,s)0 2 + S 4 -r X 2 (11) 
dx 

^ = -2r0 X (12) 
ax 

The potential V(4>,x) is given by 
1 



V(<j>,x) 



If + 



It can be projected in the x — direction to give 



1 



^,o) = -(^-i)v-ir 



(14) 



It is of the eight order power in <j> and admits symme- 
try breaking for s positive and negative. It can also be 
projected in the x — direction; in this case it gives 



V{Q,x) = \{r X 2 -lf 



(15) 



which is of the fourth order power in x an d admits sym- 
metry breaking for r positive. The model contains the 
pair of minima v± t 2 — (±1, 0) and also 113,4 = (±-^/l/s, 0) 
if s is positive, and also W5.6 = (0, ± \J 1 /r) if r is positive. 

We use the potential to calculate the mass matrix, and 
to write the masses of the two fields at the diverse minima 
of the potential. They are: m| = 4(1 — s) 2 , m 2 = 4r 2 

4(1 - ,s) 2 / S , 




4r 2 /s 



at the minima (±1,0), 

at (±y/l/a,0), and m| = 4r, to 2 = Ar at 
This model admits the presence of domain walls having 
internal structure. We examine this possibility in the case 
r > and s < 0, where we have the same four minima 
of the former model. Here, for x = 0, in the BPS sector 
defined by the minima Vx t 2 = (±1; 0) the energy is given 
by (4/3) (1 + |s|/5); for cf> = 0, in the non BPS sector 
defined by = (0, ±y/l/r) the energy is 4/3y / r. Thus, 
the ratio (5+\s\)/5y/r shows that one favors the presence 
of the BPS defect with x = as the host defect for y/r > 
l + |s|/5, and for y/r < l + |,s|/5 the host defect is the non 
BPS defect with = 0. Thus, for r G (0, 1) and for s < 
the host defect is necessaril y th e non BPS defect that 
connects the minima (0, ±y/l/r). And in this case the 
4> meson prefers to live outside the host defect, while the 
X meson prefers to live inside. This situation is similar 
to the one that appears in the former model. However, 
for r > 1 and for s < 0, if one considers y/r > 1 + |s|/5 
the host defect becomes the BPS defect that connects 
the minima (±1, 0). In this case, the x mesons prefer to 
live inside the host defect. But for the 4> mesons we now 
have two possibilities: they prefer to live outside the host 
defect if y/r > 1 + \s\, and for 1 + |s|/5 < \/r < 1 + \s\ 
they also prefer to live inside the host defect. The present 
model is more general than the former one, giving rise to 
the case where the host defect entraps both the <f> and x 
mesons. 

We notice that W in Eq. ( |Io|) is not harmonic, so the 
model is not exactly solvable. Moreover, since we have 
been unable to find an integrating factor for Eqs. ( |ll| ) 
and @, we could not be sure that the BPS states of 
the model can be solved exactly. For this reason, let us 
now use the trial orbit method to explore the presence 
of BPS states in this model. Because this new model 
contains an extra term, of the fourth order power in cf>, 
we try the orbit a + b<fi 2 + cx 2 + dcf) 4 = 0. We follow the 
former steps to get to the good orbit rx 2 = s(l — <fi 2 ) 2 , 
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for solutions connecting the minima (±1,0), with energy 
(4/3) 1 1 - s/5|. We use this orbit to rewrite Eq. Q) 
as d(p/dx = (1 — s)(l — 2 ), which is solved by cp(x) 
'■ii i — s)x}. We use this orbit to get for the other field 
the solutions 



X(%) = ± 



cosh [(1 — s)x] 



(16) 



which requires that s/r > 0. This pair of solutions exists 
for r and s positive or negative, so it appears in the sector 
connecting the minima (±1, 0) despite the presence of the 
other minima. It goes to the solution which describes a 
straight line orbit in the former model in the limit s — > 0. 

Let us now suppose that s is positive, so that there 
are minima at v^fi = (±y/l/s, 0). In this case we can 
find BPS solutions connecting those minima, with energy 
(4/3y/s)\l — l/5s|. We follow the same procedure to find 
the orbit rx 2 = (1 — s(f> 2 ) 2 / s. This orbit allows obtaining 
the solutions 



4>{x) = s/l/s tanh[^/l/s (s - l)x] (17) 

(18) 



X(x) = ± 



cosh 



'1/s (s- l)x 



which requires r > 0, so that there are minima also at 
(0, iyT/r). This last solution appears if both r and s 
are positive, so that the model must contains all the six 
minima (±1,0), (±y/TJs), and (Q,±y/\Jr). 

We further illustrate the trial orbit method investiga- 
timg another model, defined by 



1 r 

— r — 

2 6 



ls_ 

3^ 



(19) 



where r and s are real and positive, s €E (0,1). The 
presence of interactions leading to negative power in the 
fields is not unusual. They appear in several different 
contexts, for instance in vortices in planar Abelian sys- 
tems involving generalized permeabilities (2^, in ef- 
fective Yang-Mills theories coupled to scalar field, with 
appropriate color dielectric function to mimic quark and 
antiquark interactions |^5|, |26|, |2^], and in models used 
to map biological systems, involving an activator and its 
antagonist, the inhibitor pS[ ]. 

In this model the first order equations are 



dx 

dx 
dx 



1 X 
— r — 

2 6 



X 
r — 



(20) 
(21) 



We investigate the BPS sector defined by the minima 
(±1,0). Its energy is given by (4/3)(l — s). In this model, 
the presence of interactions leading to negative power 
in the <j> field makes the search for defects harder than 
before, so this example offers another good illustration 
of the trial orbit method. 

We follow the former steps and try the orbit a + bcj) 2 + 
ccj) 4 + d(/) 2 x 2 = 0. It gives the good orbit 4> 2 X 2 ~ 2s(l — 



2 ) = for r = 1. We use this orbit in the first order 
"(I) to get d<p/dx — 1 — <j> 2 , which is solved by 



equation 
4>{x) = tanh(.x 



This means that the other field x 1S 



X(x) = ± 



'2s 



sinh(x) 



(22) 



We see that xi x ) diverges for x — > 0, spliting the orbit 
into two segments. Because of this singularity we can- 
not calculate the corresponding energy, and this indicates 
that this pair of solutions is unstable. 

The last example is obtained as a two-field generaliza- 
tion of a model introduced in Ref. [^9| , which describes a 
vacuumless potential. The potential has a maximum at 
= 0, and monotonically decreases to zero at <f> — * oo. 
Potentials with such asymptotic behavior can arise in 
supersymmetric gauge theories []30[ , and more recently 
they have been used in the form of quintessence models 
fsif . They support defects which present quite different 
properties and evolution, as compared to the usual de- 
fects p9| . To generalize the model introduced in |p9| to 
the case of two real scalar fields we follow Ref. [|32|. We 
introduce the superpotential 



W{(j>, x) = arctan[sinh(^)] + W(<f>, x) (23) 



where 



WU )= /i"Msx)-sinh4>) 
v cosh(0) cosh(sx) 



(24) 



with r and s as real parameters. The first order equations 



1 



dx 

dx 
dx 



cosh(0) 



1 + sinh(^) sinh(sx) 



cosh ((/)) cosh(sx) 



1 + sinh(</>) sinh(sx) 



cosh(0) cosh (s%) 



(25) 
(26) 



These equations are much harder to solve, but we can 
use the trial orbit method to see that the orbit <f> = sx is 
a good orbit for r = 1/(1 + s 2 ). In this case we get the 
solutions 



(x) = sx(x) = arcsinh 



2 

s x 



(1 



(27) 



This pair of solutions is of direct interest to high energy 
physics, as we can see for instance in [^9[ p3| , and in Q 
in the case of quintessence with coupled scalar fields. 

The above examples illustrate how efficiently the trial 
orbit method can be used to obtain explicit BPS solutions 
in specific models, and this will certainly help exploring 
other systems, involving the yet unknown possibilities of 
finding BPS states in models described by three or more 
real scalar fields, with discrete symmetry. 
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